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ABSTRACT 

We present data from a Chandra observation of the nearby cluster of galaxies Abell 576. The core of the cluster 
shows a significant departure from dynamical equilibrium. We show that this core gas is most likely the remnant 
of a merging subcluster, which has been stripped of much of its gas, depositing a stream of gas behind it in the 
main cluster. The unstripped remnant of the subcluster is characterized by a different temperature, density and 
metalicity than that of the surrounding main cluster, suggesting its distinct origin. Continual dissipation of the 
kinetic energy of this minor merger may be sufficient to counteract most cooling in the main cluster over the 
lifetime of the merger event. 

Subject headings: cooling flows — galaxies: clusters: individual (Abell 576) — intergalactic medium — 
X-rays: galaxies: clusters 



1. INTRODUCTION 

The study of the cores of clusters of galaxies has 
undergone a renaissance in the past few years with 
the launch of the Chandra and XMM-Newton observa- 
tories. While previous observatories lacked the spatial 
resolution necessary to resolve structure within the cores 
of clusters, this new generation of telescopes has re- 
vealed an astonishing level of complexity in the struc- 
ture of the intracluster medium (ICM). Many clusters 
previously thought to be relaxed, regular systems have 
proven to be far from dynamical equilibrium, particularly 
in thei r cores (e.g.[M arkevitch. Vikhlinin, & Mazzottal l2001t 
Mazzotta. Edge. & Mar kevitchl 120031) . Abell 576 is no ex- 
ception. Earlier data, especially optical spectra of the clus- 
ter's galaxy popu lation, hinted at dynamical complexity in 
the cluster core tMohr et al.l 119961) . but earlier X-ray ob- 
servations showed the cluster to be quite regular and even 
cooler in the center, suggesting either a cooling flow or multi- 
phase gas with a v ery cool component ( Rothenflug et alJl984t 
Moh retalll996l) . 

With its low redshift, Abell 576 makes excellent use of the 
capabilities of Chandra, allowing us to examine in detail the 
very core of the cluster. In this paper, we focus on the dynam- 
ical activity in the core of cluster. The cluster shows strong 
evidence, first suggested by Mohr et al. ( 1996) from an anal- 
ysis of the galaxy population, of the remnant core of a small 
merged subcluster. We demonstrate that the X-ray data are 
consistent with this picture, and even suggest it as the most 
likely origin for the non-equilibrium gas at the center of Abell 
576. In fact, the subcluster may still be in the process of set- 
tling into the center of the main cluster's potential. 

Throughout this paper, we use the cosmologi cal parameters 
deriv ed from the first release WMAP results iBennett et alJ 
2003), so 1" = 0.738 kpc at z= 0.0377. All errors are quoted 
at 90% confidence unless otherwise stated. 

2. OBSERVATION AND DATA REDUCTION 

The data were obtained during Chandra Cycle 3 in a sin- 
gle exposure of 38.6 kiloseconds. The focus was set on the 
back-illuminated S3 chip, although significant flux from the 
cluster is detected on the adjacent front-illuminated S2 chip. 



The standard background reduction for Very Faint ( VF) mode 
data was applied using the CIAO tool acis_process_events. 
In addition, the data were corrected for charge transfer in- 
efficiency (CTI) using the CXC/MIT CTI-corrector in CIAO 
version 2.3 1 , although this only corrects the front-illuminated 
chips such as the S2. They were then filtered on the stan- 
dard ASCA grades 0, 2, 3, 4, and 6. The data were processed 
using version 2.18 of the Chandra CALDB. CTI-corrected 
blank sky background files 2 provided by Maxim Markevitch 
were used for the background correction. The VF mode back- 
ground reduction was also applied to the background files. 

The final 10.7 ksec of the observation was contaminated by 
a background flare. The flare was of the hard spectrum vari- 
ety, which affects both the front- and back-illuminated chips. 
Unlike the soft flares that Chandra sometimes experiences, 
which affect only the back-illuminated chips, the spectra of 
these hard flares are not consisten t from flare to flare, and 
therefore cannot be modeled (Markevtich 2002). We there- 
fore excised the contaminated data, and only consider here 
the first 27.9 ksec of the observation. 

Because of the small field of view of the S3 chip and the 
low redshift of Abell 576, our analysis is mostly restricted to 
the central ^200 kpc of the cluster. 

For the spectroscopic analysis we considered only data in 
the range 0.5-8 keV. Below this range the calibration of the 
ACIS CCDs is less certain, and above this range the data 
contain few photons and are dominated by particle back- 
ground. Spectral response matrices were corrected for the re- 
duction in quantum efficiency at low energies using the time- 
dependent but spatially invariant acisabs model provided by 
the Chandra X-ray Center. Because of spatial variations in 
the contaminant, the correction is not exact, especially be- 
low 1.2 keV where the effect increases. For this reason, we 
kept the absorption column fixed to the Ga lactic value of 
5.71 x 10 2() cm" 2 (iDickev & Lockmanlfl990l) in the spectral 
fits presented here. As a test, we allowed the absorption col- 
umn to vary, and measured a column density nearly identical 
to the Galactic value, but with substantial errors, typically on 
the order of 10%. The only effect of allowing the absorption 
to be a variable parameter was to increase the errors in our 
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FIG. 1. — Gaussian smoothed, exposure-corrected, 0.3-6.0 keV image of 
Abell 576. Both the S2 and S3 chips are shown. The north and southeast 
surface brightness edges are clearly visible; the west edge is less distinct. 



determinations of the other free parameters. Thus, we con- 
cluded that it was preferable to leave the absorption fixed at 
the Galactic value. 

3. BRIGHTNESS EDGES 

Figure ^ shows a Gaussian smoothed, exposure-corrected 
image of the cluster. At least two, and perhaps more, sur- 
face brightness edges are visible within the central 50 kpc. 
As shown below, they encompass a region of cool, high- 
metalicity gas. This cool gas also extends in a finger to the 
north of the cluster core, slightly west of center. As we will 
discuss in the remainder of this section, we believe this finger 
of gas to have originated in a small subcluster which is cur- 
rently accreting into the center of the main cluster. The orien- 
tation of the edges are not consistent with gas simply sloshing 
back and forth in a m ore or less fixed potential, which would 
create parallel edges (Ouilis, Bower, & Baloah 2001) as op- 
posed to the roughly triangular configuration observed. The 
observed edges are more consistent with being the outer edges 
of a wake of stripped gas left behind by a merging subcluster. 
In this picture, the subcluster initially fell in from the north, 
slightly to the west of the main cluster's center, passed the 
main cluster center once, and is now making it's second pass 
of the cluster's center. The west and southeast edges describe 
the outer edges of the wake of stripped gas from it's previous 
and current pass of the main cluster's center. This hypothe- 
sis also neatly explains the finger of gas to the north, which 
cannot be easily explained by simple sloshing. 

We extracted surface brightness and spectral profiles across 
all three edges using the regions shown in Figure |2] The 
brightness edge 40" north of the peak of the X-ray emission 
shows by far the largest jump in surface brightness: a factor 
of 1.8 ±0.15 (la) increase across the discontinuity (see Fig- 
ure[3}. A large jump in the abundance is also visible across 
the discontinuity, while the temperature does not change sig- 
nificantly (see the points with red error bars in Figure©. This 
jump in abundance is significant at more than 90% confidence 
across the north edge. At the low temperature of Abell 576, 
the increased abundance across the edge has a non-negligible 




FIG. 2. — Regions used for the spectral and surface brightness profiles. 
The exact regions used for the spectral profiles are indicated. The surface 
brightness profiles used the same regions, but subdivided into smaller radial 
bins. 




FIG. 3. — Surface brightness profile across the north edge. The solid line 
ou tside 30 kpc is a beta model fit to the data using the parameters derived 
by Mohr et al. 1 1995) using the full field of view of the Einstein IPC. Inside 
30 kpc, three models for the surface brightness jump are shown: the solid 
line indicates the same beta model after adjusting for the enhancement due 
to both the density and abundance increases inside the edge. The abundance 
only and density only components are shown as the dashed and dotted lines, 
respectively . For the /3-model, /3 = 0.64 and r c = 169 kpc. The core radius 
of the model is indicated as "r c " for reference. The errors bars are la. 



effect on the emissivity of the gas. This is illustrated in Fig- 
ure |5J which shows the surface brightness profile across the 
north edge. Outside the brightness edge, the solid line is a 
/3-model fit to the data using the core ra dius and slope de - 
termined from observations with Einstein (Mohr et al. 1996). 
Inside the edge, we use the same /3-model, but we increase the 
emissivity by an amount expected from each of three different 
models: the dotted line indicates the increased surface bright- 
ness due to the increase in density alone; the dashed line in- 
dicates the increased emissivity due to the higher abundance; 
and the solid line is the increased brightness due to both ef- 
fects. For all three models we assume spherical symmetry 
for consistency with the deprojection analysis. As the figure 
demonstrates, neither the added emissivity from the higher 
abundance nor that from the increased density can account 
entirely for the observed increase in surface brightness, but 
the two effects combined reproduce the overall normalization 
of the central brightness quite well. We note that had we kept 
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kT (Z) M kT mia kT 2 

model (keV) (Z ) (M yr" 1 ) (keV) (keV) \ 2 d.o.f. 

MEKAL 4.05^g - 3 2!no7 '■■ ■■■ ■■■ 357 297 

mekal + mkcflow 4.91™fj 0.37™|j° 5.0+|,> 0.01!°^ ■■■ 322 295 

mekal + mekal 4£7^J 0A7 ^nl 0.25+° j 289 295 
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FIG. 4. — Spectral profiles in three sectors to the north, southeast, and west 
from the cluster center. The top set of points are temperatures; the bottom set 
are abundances. All error bars are 90% confidence. 



the abundance in our spectral models fixed, we would have 
underestimated the emissivity of the gas inside the bright- 
ness edge, and would therefore have overestimated its density. 
This would then overestimate the density contrast across the 
edge, which would negatively affect the analysis that follows. 

We deprojected the surface brightness, which, when com- 
bined with the temperature and abundance profiles, allowed us 
to determine deprojected density and pressure profiles across 
the edge, under the assumption of spherical symmetry. We 
find a small jump in both the density and the pressure across 
the north edge. The density increases by a factor of 2.8^'f 
while the pressure increases by a factor of 2.4 ±0.8 (both 
lcr). In order for the higher density gas to remain confined, 
the pressure difference across the edge must be balanced by 
ram pressure from motion of the high density gas through the 
lower density gas. The observed pressure difference implies 
that the higher density gas is moving through the lower den- 
sity gas with a velocity of 750 ±270 km s" 1 , or Mach0.9±0.3 
at the sound speed of the lower density gas (both errors are 
lcr). We note that if we had failed to account for the increase 
in abundance across the edge, that is, if we had assumed a 
constant abundance on both sides of the edge, we would have 
overestimated the density jump across the edge and conse- 
quently would have overestimated the velocity of the dense 
gas cloud. The velocity we measure is consistent with veloc- 
ities of both merging/accreting subclusters measured in other 
clusters (e.g. Markevitch et al. 2000) and with velocities mea- 
sured for some "slo shing" edges in the cores of relaxed cluster 
( Markevtich 2003). On its own, then, the measured velocity 
of the north edge is incapable of distinguishing between these 
two scenarios for the creation of the non-hydrostatic features 
in the cluster core. 

To the southeast of the cluster center, a fainter edge is vis- 
ible in the image (see Figure Another yet fainter edge 



appears to the west. Both of these edges display the same 
abundance gradient as the north edge, though in both cases 
the abundance jump appears to be more of a gradient than a 
sharp edge, and is measured with much less significance (see 
Figure^. 

A combined spectral analysis of all 3 of sectors yields im- 
proved statistics at the cost of spatial resolution. This is plot- 
ted as the solid grey boxes in Figure The best-fit models 
for the combined spectra outside the edges are given in Ta- 
ble ^ The columns are as follows: (1) spectral model; (2) 
& (3) temperature and abundance, respectively, of primary 
MEKAL model; (4) Cooling rate of MKCFLOW model; (5) 
minimum temperature of cool gas in MKCFLOW model; (6) 
temperature of second MEKAL model, where the abundance 
of the second model is tied to that of the primary MEKAL 
model; (7) & (8) \ 2 ar, d number of degrees of freedom for 
the fit. From this combined analysis, we find that the central 
abundance is different from the abundance outside the edges 
at greater than 90% confidence. As expected, the tempera- 
ture measured outside the edges is approximately equal to the 
mean temperature measured with higher spatial resolution in 
the individual sectors. However, when fit with a two tem- 
perature model, the best-fit high temperature is only slightly 
higher than for the single temperature model, while the cool 
temperature is 0.25^ \ keV and contributes only 5% to nor- 
malization of the model spectrum. (Because of the errors on 
this measurement and given that our spectral fit only includes 
data at energies > 0.5 keV, this cool component should be 
considered an upper limit of 0.65 keV.) Similarly, a single 
temperature plus cooling flow model yields a relatively small 
mass accretion rate of 5 ± 1 M Q yr" 1 but to a very low tem- 
perature of less than 0.2 keV. An F-test shows that the two 
temperature model is a better fit than the single temperature 
model at about the 3er level, while the cooling flow model is 
better than the single temperature by only ~ 1 a. In both multi- 
phase models, the cooler component contributes a very small 
fraction of the total emission. 

While a radial temperature gradient across the region used 
for our analysis (~30-150kpc) would show qualitatively sim- 
ilar effects, we would expect such a gradient to cover a rela- 
tively small range of temperatures, say, a facto r of 2 as is seen 
in the most relaxed cooling flow systems (e. g. iPeterson et alJ 
l200lHBlanton. Sarazin. & McNamaral2003l) . To find gas at a 
temperature less than 1/4 of the ambient temperature is quite 
unusual. The merging subcluster hypothesis mentioned above 
provides an attractive explanation for this extremely cool gas, 
if the cool gas has been stripped from the subcluster during 
its infall. The small radius of curvature of the north edge sug- 
gests that this remnant core is physically quite small and that 
the original subcluster was also relatively small and therefore 
cool. This provides a consistent qualitative explanation for 
the presence of the cool component. Unfortunately, given the 
large amount of gas likely to have been stripped from the sub- 
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FIG. 5. — (a) Heating rate and luminosity of radiative cooling as a function of radius. The exact definition of the heating rate is given in the text. The dotted 
curves indicate the minimum and maximum allowed rates given the errors on our measurement of the subcluster's velocity, (b) Integrated cooling rate as a 
function of radius. The solid line is the rate in the absence of heating. The dotted lines are the resulting integrated cooling rate if the kinetic energy of the cool 
core is dissipated over the given timescales. The dissipated energy is assumed to heat the gas at a constant rate per volume over the entire region. 



Table 2. Spectral fit parameters for cool core 
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(keV) 
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(M Q yr" 1 ) 
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X 2 d.o.f. 



mekal 3 - 48 -Qjj °- 56 ^)'n "' "' "' 186 188 

mekal + mkcflow 3.77;5il 0.64™ : j| 1.1 ±0.5 O.Ol^g; ■•• 183 186 
MEKAL + MEKAL 3 - 62 -0 1Q a66 ^)'« "' "' °- 65 -045 1 81 186 

MEKAL + FOREGROUND 3.21™^ O^ojjj'Tq 'J_ 4.05" " 186 188 

"Temperature and abundance of foreground component taken from single-temperature fit in TablefTI 



cluster during infall, it is not possible for us to determine the 
subcluster's original mass and virial temperature and is there- 
fore not possible to conclude whether or not the cooler gas 
is consistent with having originated in the undisturbed sub- 
cluster. We note, however, that the "sloshing" picture for the 
brightness edges cannot explain the existence of the cool com- 
ponent. Sloshing wo uld create the same gen eral appearance 
of a cold front (e.g. Markevitch et al. 2001), but would not 
cause the higher density gas to cool any more than it would in 
the absense of the sloshing. 

An alternative explanation for the presence of gas at 
less than 1/4 of the ambient temperature is that the usual 
method of preventing a classical cooling flow from develop- 
ing, whether that be AGN heating, conduction, or some other 
mechanism, is absent or suppressed in A576. Indeed the ab- 
senc e of the characteristic "bubbles" seen in other clusters 
(e.g. Fabia n et alJ2000HMcNamara et alJ2000llBlanton et al.1 
2001) suggests there has been no recent strong AGN activity, 
for one. However, we find this explanation to be unsatisfac- 
tory, since it would require the existence of a strong cooling 
flow, which is explicitly ruled out by the small cooling rate 
determined above. 

4. COOLING FLOW 

As Figure 0] shows, the temperature drops in the very core 
of the cluster, that is, inside the brightness edge. It is natural to 
ask, then, if this gas shows any evidence of being multi -phase. 
To test this, we fit a spectrum of the gas in the core with a 
single-temperature absorbed MEKAL model (Kaastra 1992), 
with a MEKAL model plus a multi-phase MEKAL model 
(MKCFLOW model), and with the sum of two MEKAL mod- 
els. The salient parameters of the fits are shown in Table [2] 



The columns are as follows: (1) spectral model; (2) & (3) 
temperature and abundance, respectively, of primary MEKAL 
model; (4) Cooling rate of MKCFLOW model; (5) minimum 
temperature of cool gas in MKCFLOW model; (6) tempera- 
ture of second MEKAL model, where the abundance of the 
second model is tied to that of the primary MEKAL model; 
(7) & (8) x 2 an d number of degrees of freedom for the fit. As 
shown by the similar fit statistics, the goodness of fit of all 
three models are essentially identical. 

The spectroscopic cooling rate we measure is quite low: an 
order of magnitude smaller than the "classical cooling flow" 
accretion rate of M = M/t mo \ = 11 M Q yr" 1 derived from the 
gas mass and cooling time in the inner 30 kpc (see Figure [5] 
and |5]below). We note that the cool component in the two- 
temperature model only contributes ^0. 1 % to the overall nor- 
malization of the model. This is consistent with the extremely 
small cooling rate measured for the cooling flow model. 

When the second temperature component of the two- 
MEKAL model was set to the parameters of the fit to the outer 
gas (the equivalent of deprojection), the fit was no better than 
for the single-temperature model. This is perhaps not too sur- 
prising, however, since the contribution of the projected fore- 
ground component to the overall normalization of the model 
is only about 2%. 

Surprisingly, the best-fit low temperature in the cooling 
flow model, that is, the temperature to which the gas cools, 
is equal to the minimum temperature allowed by the model, 
£7] ow = 0.01 keV. However, the error on the fit allows for kT\ ovl 
to be as high as 0.9 keV. This is about 1/4 of the ambient 
temperature, which is consistent with minim um temperature s 
found in stronger cooling flows (e.g. iPeterson et all 2001). 
The low temperature of this gas is consistent with that of the 
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cool gas found at larger radii. 

In short, we find essentially no evidence of multi-phase gas 
in the central cool core. Of the models described above, the 
cooling flow model contributes the largest contribution from 
multi -phase gas to the overall emission, with only 5% of the 
emission coming from the multi-phase gas. 

5. SUPPRESSION OF COOLING 

Given that we observe essentially no cooling in either the 
core or outside the core in what is otherwise a quite relaxed 
cluster, we now explore whether the dissipation of the kinetic 
energy of the remnant core is capable of suppressing cooling 
at the observed level. In the rest of this section we take the 
merging subcluster hypothesis to be the correct explanation 
of the brightness edges. 

We determined the kinetic energy of the gas inside the north 
edge using the velocity of the edge we measured above plus a 
gas mass determined from the deprojected density profile in- 
side the edge (~2x 10 10 M Q ). We then calculated the rate of 
energy input from the dissipation of this kinetic energy over 
a variety of timescales. Figure^ shows this energy dissipa- 
tion rate compared to the luminosity due to radiative cooling 
as a function of radius. The timescale used for calculating the 
heating rate is three crossing times of the cluster to the given 
radius at the current velocity of the north edge. In three cross- 
ing times, the moving cool core will have swept up its mass in 
gas, reducing its kinetic energy by 3 /4. We therefore assume 
perfectly efficient thermalization of 3/4 of the kinetic energy 
over a timescale equal to three crossing times by the core at 
its current velocity in calculating the heating rate. 

If we take the point at which the west and southeast edges 
converge as indicative of the current orbital radius of the sub- 
cluster (^100 kpc from the cluster center), the dissipation 
timescale derived using the above method is 4 x 10 s years. 
As can be seen from Figure [5J), the heating simply from the 
dissipation of the kinetic energy of the subcluster is capable 
of suppressing cooling by a factor of more than 4 in the inner 
100 kpc over this timescale. If the dissipation of the core's 
kinetic energy is spread out over 10 9 years, cooling can still 
be suppressed at the level observed in the inner 30 kpc. Nu- 
merical simulations have shown that mergers of unequal mass 
clusters thermalize their k inetic energy on ti mescales of one 
to a few times 10 9 years (Ricker & Sarazin 2001), although 
the mass ratio involved in A576 is probably much larger than 
has generally been tested in these simulations. 

So far, we have omitted any discussion of the role of tur- 
bulence in suppression of cooling. Turbulence generated 
by the motion of such a small subcluster would eventu- 
ally dissipate and be thermalized on a timescale < 10 8 yr 
(Fuiita, Takizawa, & Sarazin 2003), which is short compared 
to the timescale we are considering for dissipation of the ki- 
netic energy. We therefore consider it reasonable to assume 
close to 100% thermalization efficiency. We have also ne- 
glected the influence of the dark matter on the total energy 
budget. This is not so easily dismissed, as the dark matter 
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halo of the subcluster will take longer to dissipate its kinetic 
energy than will the gas. The dark matter oscillations about 
the center of mass of the system will induce additional motion 
in the gas, transferring kinetic energy from the dark matter to 
the gas and thereby increasing the total energy available for 
heating the gas. This will only serve to amplify the effect we 
have discussed. Our naively calculated heating rate should 
therefore be taken as a lower limit to the heating resulting 
from the subcluster merger. 

Other alternatives for balancing cooling in clusters have 
been discussed. These include heating from central AGN 
iChurazov et al]l2002l iBriiggen & K aiser 2002; Fabia n et"al] 
2003), from conduction ( e.g. iNa ravan &Medvedev 2001), 
and from turbulent mixing ( Kim & Naravan 2003 ) induced by 
AGN activity. Our data show no evidence for the X-ray cav- 
ities typical of clusters with strong central AGN, and indeed 
the cluster has one extremely weak AGN at its center with no 
non-nuclear radio emission. Since these cavities tend to per- 
sist fo r several AGN duty cycles (e.g. Perseus, Fabi an et"aT1 
2000) and none are observed in A576, AGN heating can be 
ruled out, at least over the last few times 10 8 years. 

6. CONCLUSIONS 

While earlier X-ray observations of Abell 576 have shown 
it to be quite regular on large scales, we have demonstrated 
that the core of the cluster is far from dynamical equilibrium. 
We found multiple surface brightness edges in the cluster cen- 
ter which we have demonstrated to be indicative of jumps in 
both density and abundance. Of the two most likely explana- 
tions for the existence of these edges, our analysis favors the 
hypothesis that they are formed by gas stripped from a merg- 
ing subcluster. Most of the gas appears to have been stripped 
from the subcluster, leaving a core only ^30 kpc in radius. 
The stripped gas has been found to have both a lower temper- 
ature and a higher abundance than the gas in the rest of Abell 
576. 

We find no evidence of gas cooling from the ambient tem- 
perature of the main cluster, but do find some suggestion of 
very cool gas at a temperature expected of gas that had con- 
densed out of the ICM of the subcluster. The simple cooling 
rate derived from the gas mass and cooling time is an order 
of magnitude larger than the spectroscopically measured M. 
We have demonstrated that dissipation of the kinetic energy 
of the observed remnant core of an infalling subgroup may be 
sufficient to reduce cooling to the observed rate, if that energy 
is dissipated over a timescale of < 10 9 years. 
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